Measurement of Angular Distributions of Drell-Yan Dimuons in p + p Interactions at 

800 GeV/c 
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We report a measurement of the angular distributions of Drell-Yan dimuons produced using an 800 
GeV/c proton beam on a hydrogen target. The polar and azimuthal angular distribution parameters 
have been extracted over the kinematic range 4.5 < < 15 GeV/c 2 (excluding the T resonance 
region), < pr < 4 GeV/c, and < if < 0.8. The p + p angular distributions are similar to those 
of p + d, and both data sets are compared with models which attribute the cos 2<f> distribution either 
to the presence of the transverse-momentum-dependent Boer-Mulders structure function hi or to 
QCD effects. The data indicate the need to include QCD effects before reliable information on the 
Boer-Mulders function can be extracted. The validity of the Lam- Tung relation in p + p Drell-Yan 
is also tested. 

PACS numbers: 13.85.Qk, 14.20.Dh, 24.85. +p, 13.88.+e 



The study of the transverse momentum dependent 
(TMD) parton distribution functions of the nucleon has 
received much attention in recent years as it provides 
new perspectives on the hadron structure and QCD 1]. 
One of these TMD distribution functions, first consid- 
ered by Sivers Q, represents the correlation between 
the quark's transverse momentum, fcj_, and the trans- 
verse spin of the nucleon, S±. This so-called Sivers func- 
tion, fij'(x ) fcjj, where x is the fraction of proton's mo- 
mentum carried by the quark, is time-reversal odd (T- 
odd) and can arise from initial- or final-state interac- 
tions [3|. More generally, the requirement of gauge in- 
variance of parton distributions was shown to provide 
nontrivial phases leading to the existence of T-odd distri- 
bution functions Recent measurements of the semi- 
inclusive deep-inelastic scattering (SIDIS) by the HER- 
MES [6] and COMPASS collaborations have shown 
clear evidence for the presence of the T-odd Sivers func- 
tions. These data also allow the first determination |8| 
of the magnitude and flavor structure of the Sivers func- 
tions. 



Another T-odd distribution function is the Boer- 
Mulders function, hi(x, kj_), which signifies the corre- 
lation between k± and the quark transverse spin, s±, in 
an unpolarized nucleon . The Boer-Mulders function is 
the chiral-odd analog of the Sivers function and also owes 
its existence to the presence of initial/final state interac- 
tions While the Sivers function is beginning to be 
quantitatively determined from the SIDIS experiments, 
very little is known about the Boer-Mulders function so 
far. 

Several model calculations have been carried out for 
the Boer-Mulders functions. In the quark-diquark model, 
it was shown that the Boer-Mulders functions are identi- 
cal to the Sivers functions when only the scalar diquark 
configuration is considered [llT |. More recently, cal- 
culations taking into account both the scalar and the 
axial-vector diquark configurations found significant dif- 
ferences in flavor dependence between the Sivers and 
Boer-Mulders functions In particular, the u and 

d valence quark Boer-Mulders functions are predicted to 
be both negative, while the Sivers function is negative 



2 



for the u and positive for the d valence quarks. Other 
calculations using the MIT bag model [l]| , the relativis- 
ts constituent quark model 14 1 , the large- N c model [l5| , 
and lattice QCD [l6| also predict negative signs for the 
u and d valence Boer-Mulders functions. Burkardt re- 
cently pointed out [ijj that the negative signs for the 
Boer-Mulders functions are expected for both nucleons 
and pions. The model predictions for the same signs of 
the u and d Boer-Mulders functions remain to be tested 
experimentally. Furthermore, the striking prediction Q 
that the T-odd Boer-Mulders functions in the SIDIS pro- 
cess will change their signs for the Drell-Yan process also 
awaits experimental confirmation. 

The Boer-Mulders functions can be extracted from 
the azimuthal angular distributions in the unpolarized 
Drell-Yan process, \i\h2 — * l + l~x. The general expres- 
sion for the Drell-Yan angular distribution is [19( 



da 7 v „ 

-j^j oc 1 + A cos 6 + (i sin 26 cos <p + — sin 



'cos 20, (1) 



where 9 and <ft are the polar and azimuthal decay an- 
gle of the l + in the dilepton rest frame. Boer showed 
that the cos 20 term is proportional to the convolution of 
the quark and antiquark Boer-Mulders functions in the 
projectile and target This can be understood by 

noting that the Drell-Yan cross section depends on the 
transverse spins of the annihilating quark and antiquark. 
Therefore, a correlation between the transverse spin and 
the transverse momentum of the quark, as represented 
by the Bocr-Muldcrs function, would lead to a preferred 
transverse momentum direction. 

Pronounced cos 20 dependences were indeed observed 
in the NA10 [H and E615 21] pion-induced Drell-Yan 
experiments, and attributed to the Boer-Mulders func- 
tion. The first measurement of the cos 20 dependence 
of the proton-induced Drell-Yan process was recently re- 
ported for p + d interactions at 800 GeV/c [22l |. In 
contrast to pion-induced Drell-Yan, significantly smaller 
(but non-zero) cos20 azimuthal angular dependence was 
observed in the p + d reaction. While the pion-induced 
Drell-Yan process is dominated by annihilation between a 
valence antiquark in the pion and a valence quark in the 
nucleon, the proton-induced Drell-Yan process involves 
a valence quark in the proton annihilating with a sea 
antiquark in the nucleon. Therefore, the p + d result 
suggests [22| that the Boer-Mulders functions for sea an- 
tiquarks are significantly smaller than those for valence 
quarks. 

A recent analysis 23] indicated that the E866 p + d 
data are consistent with the u and d Boer-Mulders func- 
tions having the same signs, as predicted by various mod- 
els. However, the p + d data alone cannot provide an 
unambiguous determination of the flavor dependence of 
the Boer-Mulders functions. Moreover, it was recently 
pointed out 24 , [2|| that QCD processes would lead to 
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FIG. 1: Parameters \,[i,v and 2p — (1 — A) vs. pr in the 
Collins-Soper frame. Solid squares (open circles) are for E866 
p + p (p + d) at 800 GeV/c. The vertical error bars include 
the statistical uncertainties only. 



account in the extractions [3, [2^, [26[ of Boer-Mulders 
functions from the Drell-Yan data. In this paper we re- 
port the Drell-Yan angular distributions of the p+p reac- 
tion at 800 GeV/c, which provides further constraints on 
the flavor dependence of the Boer-Mulders functions [2(| . 
We also compare the cos 2(f) dependences of p+p and p+d 
data with the prediction of QCD. 

The Fermilab E866 experiment was performed us- 
ing the upgraded Meson-East magnetic pair spectrom- 
eter 



27]. An 800 GeV/c primary proton beam with up 



a sizeable cos 20 effect which has not been taken into 



to 2 x 10 12 protons per 20 s beam spill was incident upon 
one of three identical 50.8 cm long target flasks contain- 
ing either liquid hydrogen, liquid deuterium or vacuum. 
A copper beam dump located inside the second dipole 
magnet (SM12) absorbed protons that passed through 
the target. Downstream of the beam dump was an ab- 
sorber wall that removed hadrons produced in the target 
and the beam dump. 

Several settings of the currents in the three dipole mag- 
nets (SM0, SM12, SM3) were used in order to optimize 
acceptance for different dimuon mass regions. Data col- 
lected with the "low mass" and "high mass" settings 27] 
on liquid hydrogen and empty targets were used in this 
analysis. The detector system consisted of four track- 
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TABLE I: Mean values of the A, v parameters and the quan- 
tity 2v — (1 — A) for the p + p, p + d, and ir~ + W Drell-Yan 
measurements. 





p + p 
800 GcV/c 
(E866) 


P + d 
800 GeV/c 
(E866) 


n- + W 
194 GeV/c 
(NA10) 


(A) 


0.85 ±0.10 


1.07 ±0.07 


0.83 ±0.04 


<^> 


-0.026 ±0.019 


0.003 ±0.013 


0.008 ±0.010 




0.040 ± 0.015 


0.027 ±0.010 


0.091 ± 0.009 


(21/ -(1- A)) 


-0.07 ±0.10 


0.12 ±0.07 


0.01 ±0.04 



ing stations and a momentum analyzing magnet (SM3). 
Tracks reconstructed by the drift chambers were extrapo- 
lated to the target using the momentum determined from 
the bend angle in SM3. The target position was used to 
refine the parameters of each muon track. 

From the momenta of the p + and , kinematic vari- 
ables of the dimuons (iF,m w , and pt, where xf is the 
fraction of the cm. momentum carried by dimuon of 
mass Tn^fj,, and px is the dimuon transverse momen- 
tum) were readily reconstructed. The muon angles 9 
and 4> in the Collins-Soper frame [28| were also calcu- 
lated. To eliminate the J/\& and T resonance back- 
ground, dimuon events with m w < 4.5 GeV/c 2 and 9.0 
GeV/c 2 < m MM < 10.7 GeV/c 2 were rejected in the anal- 
ysis. A total of «54,000 p + p Drell-Yan events cov- 
ering the decay angular range —0.5 < cos 9 < 0.5 and 
— 7r < <f> < ir remain. Detailed Monte Carlo simulations 
of the experiment using the MRST98 parton distribu- 
tion functions [29j for NLO Drell-Yan cross sections have 
shown good agreement with the data for a variety of mea- 
sured quantities. 

Figure 1 shows the angular distribution parameters 
A,/i, and v vs. pt- To extract these parameters, the 
Drell-Yan data were grouped into 5 bins in cos (9 and 8 
bins in <p for each px bin. A least-squares fit to the data 
using Eq. 1 to describe the angular distribution was per- 
formed. The extracted values of A, /i, v are insensitive to 
their values used in the Monte Carlo simulation. Only 
statistical errors are shown in Fig. 1. The primary contri- 
butions to the systematic errors are the uncertainties of 
the incident beam angles on target. Analysis performed 
by allowing the beam angles to vary within their ranges 
of uncertainty has shown that the systematic errors are 
small compared to the statistical errors. The E866 p + d 
Drell-Yan data are also shown in Fig. 1 for comparison 
with the E866 p + p data. The p + d data contain a total 
of «118,000 events covering an identical cos 9 range. The 
PT-averaged values of (A), (p), and (v) for p + p, p + d, 
and the NA10 n~ + W data ^ are listed in Table I. 
Within statistics, the angular distributions of p + p are 
consistent with those of p + d. Also shown in Fig. 1 and 



Table I is the quantity 2v — (1 — A), which should vanish 
if the Lam- Tung relation is valid. While QCD effects can 
lead to A ^ 1 and fx, v ^ 0, Lam and Tung showed [IB] 
that the relation 1 — A = 2v is largely unaffected by QCD 
corrections. Table I shows that while (A) deviates from 
1 and (v) is nonzero for the E866 p+p and the NA10 
7r~ ± W Drell-Yan data, the Lam- Tung relation is indeed 
quite well satisfied within statistical uncertainty for all 
Pt- This differs from the observation of a significant vio- 
lation of the Lam- Tung relation at large px by the E615 
collaboration in the 7r" ± W reaction at 252 GeV/c [2lT |. 

Figure 2 shows the parameter v vs. pt for the p + p 
and p + d Drell-Yan data. The solid curves are calcula- 
tions 23|, 2(| for p + p and p + d using parametrizations 
of the Boer-Mulders functions deduced from a fit to the 
p + d Drell-Yan data. The predicted larger values of v 
for p+p compared to p + d in the region of pt ~ 1.5 
GeV/c are not observed (the predicted p + p/p + d ra- 
tio, R, for 0.5 < p T < 2.0 GeV/c, is ~ 2, while the 
data give R = 1.0 ± 0.5). Furthermore, the shape of the 
predicted pt dependence differs from that of the data, 
resulting in a reduced x 2 value of 3.2 for 5 degrees of 
freedom (probability of 0.7%). This strongly suggests 
that there could be other mechanisms contributing to the 
cos 20 azimuthal angular dependence at large pt- In re- 
cent papers [3,[23|, the QCD contribution to the cos 2<p 
azimuthal angular dependence is given as 



JQ 2 



1 



Ql/Q 2 



(2) 



where Q± is the dimuon transverse momentum. The pre- 
dicted QCD contribution, the same for p+p and p+d due 
to the identical kinematic coverage for the two reactions, 
is shown as the dot-dashed curve in Fig. 2. A compar- 
ison between the QCD prediction with the data gives a 
reduced \ 2 °f 1-0 f° r 5 degrees of freedom (probability 
of 42%) for p + p and a reduced \ 2 °f 1-9 (probability of 
9%) for p + d. From Fig. 2 it is evident that the QCD 
contribution is expected to become more important at 
high pt while the Boer-Mulders functions contribute pri- 
marily at lower pt- An analysis combining both effects 
is required in order to extract reliably the Boer-Mulders 
functions from the p+p and p+d data. It is worth noting 
that the tt~ + W Drell-Yan data [2(| EOJ also show large 
values of v at large pt, consistent with the presence of 
QCD effects. 

The p + p Drell-Yan angular distributions have also 
been analyzed for other kinematic variables. Figure 3 
shows the values of v vs. m^^Fjii, and X2, where x\ 
and X2 are the Bjorken-a; for the beam and target par- 
tons, respectively. Again, for each bin the data were 
divided into 5 bins in cos 9 and 8 bins in <j) in order to ex- 
tract the angular distribution parameters. The p+d data 
are also shown for comparison. Figure 3 shows that the 
magnitude of v for p + p is consistent with that for p + d 
for most of the kinematic regimes. These data provide 
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is not observed in the p + p reaction. The Lam- Tung 
relation remains valid for the p + p Drell-Yan data. The 
overall magnitude of the cos 2<j> dependence for p + p is 
consistent with, but slightly larger than that of p + d. 
The data suggest the presence of higher-order QCD cor- 
rections at high px, and it is important to take this con- 
tribution into account before reliable extraction of the 
Boer-Mulders functions could be obtained. 
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gelsang, and Jianwei Qiu. This work was supported in 
part by the U.S. Department of Energy and the National 
Science Foundation. 



FIG. 2: (color online). Parameter v vs. pr in the Collins- 
Soper frame for the p + p and p + d Drell-Yan data. The solid 
and dotted curves are calculations [23I ] for p + p and p + d, 
respectively, using parametrizations based on a fit to the p + d 
data. The dot-dashed curve is the contribution from the QCD 
process (Eq. 2). 
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FIG. 3: Parameter v vs. m M(J , xf, xi, and X2 in the Collins- 
Soper frame for p + p (solid squares) and p + d (open circles) 
at 800 GeV/c. The vertical error bars correspond to the sta- 
tistical uncertainties only. 



further input for future extraction of the Boer-Mulders 
functions. 

In summary, we report a measurement of the angu- 
lar distributions of Drell-Yan dimuons for p + p at 800 
GeV/c. The pronounced cos 2<p azimuthal angular de- 
pendence observed previously in pion-induced Drell-Yan 
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